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Abstract 

Some time ago Salmhofer demonstrated the equivalence of the strong cou- 
pling lattice Schwinger model with Wilson fermions to a certain 8- vertex model 
which can be understood as a self-avoiding loop model on the square lattice 
with bending rigidity 77 = 1/2 and monomer weight z = (2k)~ 2 . The present 
paper applies two approximate analytical methods to the investigation of crit- 
ical properties of the self-avoiding loop model with variable bending rigidity, 
discusses their validity and makes comparison with known MC results. One 
method is based on the independent loop approximation used in the literature 
for studying phase transitions in polymers, liquid helium and cosmic strings. 
The second method relies on the known exact solution of the self-avoiding loop 
model with bending rigidity rj = l/\/2. The present investigation confirms re- 
cent findings that the strong coupling lattice Schwinger model becomes critical 
for K cr ~ 0.38 — 0.39. The phase transition is of second order and lies in the 
Ising model universality class. Finally, the central charge of the strong cou- 
pling Schwinger model at criticality is discussed and predicted to be c = 1/2. 
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1 Introduction 



Recently, the strong coupling lattice Schwinger model with Wilson fermions (Nf = 
1) has received some attention 0-fl following work by Salmhofer || who has shown 
that it is equivalent to a certain 8-vertex model (a 7-vertex model, more precisely) 
which also can be understood as a self-avoiding loop model on the square lattice 
with bending rigidity 77 = 1/2 and monomer weight z = (2k)~ 2 . Beyond its toy 
character interest in the lattice Schwinger model (QED 2 ) mainly derives from the 
similarity of some of its major features with those of QCD in 4D. However, because 
the result of Salmhofer |J is related to the polymer (hopping parameter) expansion 
of the fermion determinant M, || the strong coupling Schwinger model is also in- 
teresting from the point of view of the dynamical fermion problem within lattice 
gauge theory. While some investigations have been devoted to the polymer expan- 
sion of the fermion determinant in the case of staggered fermions |[L3|| to date 
almost no attention has been paid to the corresponding case of Wilson fermions 
TM due to the additional difficulties involved in general. However, while in the 



strong coupling limit the system with staggered fermions (QCD, QED) reduces to 



a pure monomer-dimer system W5 the same is not true for Wilson fermions as the 
investigation of Salmhofer || demonstrates. The equivalence of the strong coupling 
lattice Schwinger model with Wilson fermions to a self-avoiding loop model enables 
certain methods used in other branches of physics, e.g. in condensed matter physics 
(polymers, defect mediated phase transitions) and in cosmic string physics, to be 
exploited in its investigation ||16||-||21|| ■ At the same time, its equivalence (in another 
language) to some 8-vertex model [2^] makes further results available. 



Self-avoiding loop models [g3| , |24j have a long history due to their prominent role 
in polymer physics as well as their attraction as a simple problem of non-Markovian 
nature. In addition, systems of closed non-crossing lines or systems which can be 
approximated by them appear in a variety of contexts ranging from condensed mat- 
ter physics through cosmology to quantum field theory which generates common 
interest for appropriate model building PSQ , ]T5| . Recently, quantum field theoretic 
methods have been exploited to study the critical behaviour of self-avoiding loop 
models in two dimensions |26| ||29|| . Somewhat lesser attention has been paid so far 
to the self-avoiding loop model with variable bending rigidity (while for open chains 
with bending rigidity a number of investigations exists, e.g., [|30] and references 
therein). Beyond the work of Muser and Rys [|3T| certain insight in this direction 
has been obtained in connection with the study of 2D vesicles [p2| - p|| . 

From the point of view of the 8-vertex model, a general solution to the self- 
avoiding loop model with variable bending rigidity on the square lattice is not known. 
However, for the special case 77 = l/y/2 the free-fermion condition [35], [|36|1 is ful- 
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filled and it can be solved exactly [E7|]— This way one point on the critical line 
of the the self-avoiding loop model with variable bending rigidity is known exactly 
and consequently one may use methods of perturbative nature in its neighbourhood 
to approximately find the critical line there by analytical methods. 

The plan of the paper is as follows. In section 2 we shortly review the relevant 
facts concerning the lattice Schwinger model with Wilson fermions and discuss the 
relation of recent exact studies of its partition function on finite lattices |4[], || to 
the earlier MC results of Miiser and Rys [|3l|. Section 3 is devoted to the approxi- 



mate analytical study of the self-avoiding loop model (SALM) with variable bending 
rigidity by means of the independent loop approximation. Section 4 then explores 
the application of the exact solution of the SALM with bending rigidity 77 = l/\/2 
to the study of the critical behaviour of the SALM in a neigbourhood of it in the 
relevant parameter space. Section 5 finally discusses the picture emerging from the 
present investigation paying special attention to the central charge of the SALM 
with variable bending rigidity on the critical line. 



2 The strong coupling Schwinger model with 
Wilson fermions 

The partition function Za of the Schwinger model with Wilson fermions (with Wilson 
parameter r = 1) on a certain lattice A is given by the standard expression 



Z A = J DUDi/jD^ e~ b , (1) 
where D denotes the multiple integration on the lattice. The action 5* is defined by 

S = S F + PS G , (2) 
xeA V m 

+ $(x)(l - j„)Ul(x)^(x + £)) - M$(x)i/>(x)J , (3) 

and Up = exp [— iA^], M = 2 + m, (3 = 1/g 2 . So is the standard Wilson action and 
the hopping parameter k is given by k = 1/2M. Salmhofer has shown that in 
the strong (infinite) coupling limit (3 = the partition function Z\ equals that of an 
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8- vertex model (a 7- vertex model due to eq. @, more precisely) |22|] with weights 
(cf. fig. 1) 



W\ = z 



w 2 = 



4^ = M ' 



w 3 = w 4 = 1 , 

^5 = W 6 = W 7 = W 8 = 7] 

Consequently, one can write 



1 

2 



(4) 
(5) 
(6) 

(7) 



E 

L 



1 

2. 



,[A| 



1^1 „Cf(X) 



(8) 



(9) 



where L denotes any self-avoiding loop configuration, \L\ and C(L) are the number 
of links and corners respectively a polymer configuration L is built of, and |A| is 
the number of lattice points of the lattice A. Z\[z,i]\ is the partition function of a 
self-avoiding loop model (SALM) with monomer weight z and bending rigidity 77. 
The same expression can of course also be obtained for non-compact QED 2 . 

From the point of view of lattice field theory it is interesting to know the phase 
structure of the lattice Schwinger model. For free fermions {(3 = 00) the critical 
value of the hopping parameter reads K cr (/3 — 00) = 1/4. In order to pin down the 
critical line for (3 < 00 it is of particular interest to know where it ends (J3 = 0). 
There is a critical point for K cr (f3 = 0) = 00 because then the strong coupling 
Schwinger model reduces to a 6-vertex model whose behaviour is known from its 
exact solution 0, p2| . This point however is believed to be isolated and not to 
be the end point of the critical line starting at ) = 1/4 [^]. Recently, exact 

studies of the partition function of the strong coupling Schwinger model have been 
made on finite lattices [ffl, |J. It has been found K cr (0) ~ 0.38 — 0.39 and that the 
phase transition is likely a continuous one (second order or higher) |J. 

It is worthwhile to compare the result obtained in || with the MC investi- 
gation of the SALM with variable bending rigidity undertaken by Miiser and Rys 
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|31|j (see also ||t]] for some computational background). Their investigation has been 



inspired by the generalized loop model of Rys and Helfrich \24\. Muser and Rys [ [31 



employ a different parameter set than {z,r]} which we are going to describe first. 
Their language is thermodynamic in spirit and their parameters temperature and 
line stiffness {T, s} are introduced the following way 



z = e 



77 = e 



which in turn entails 



[l-s)/T 



-s/T 



(10) 
(11) 



T 



ln- 

V 

In i] 

hT^ 



(12) 
(13) 



For positive temperatures T negative values of the line stiffness s correspond to 
values of the bending rigidity 77 > 1 (i.e., bending preferred) and positive values of s 
to 7] < 1 (i.e., bending is costly). The Jacobian T of this coordinate transformation 
from {z, rj} G {[0, 00), [0, 00)} to {T, s} G {[0, 00), (—00, 00)} reads 



T 



zr\ 



T 3 



In 



s)/T 



(14) 
(15) 



and has consequently singular lines beyond the given range of the map. Figure 2 
displays the result of Muser and Rys |3T| (adapted from their fig. 2) for the critical 
line of the loop model with line stiffness and fig. 3 displays the same information in 
{z, rj} coordinates (for further comments see the figure captions). In regions II and 
III the system at criticality is found to exhibit Ising-like behaviour while in region I 
some non-universal behaviour is seen. For better orientation the ordinary loop gas 
(rj = 1) result is specially shown in figs. 2, 3 |41||-|4^|. One immediately recognizes 
that the results found for the strong coupling Schwinger model j3J , fits well onto 
the critical line given by Muser and Rys. Moreover, the MC result of Muser and 
Rys also well agrees with the exactly known critical point for the free fermion model 
(77 = 1/ V2) to be discussed in section 4. For the ordinary loop gas it has also been 
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found numerically that critical exponents (and amplitudes, in part) agree well 
with those of the Ising model. The free fermion model {rj = l/y/2, see section 4) 
of course also lies in the Ising universality class. This immediately suggest (for a 
further discussion see section 5) that in general the SALM with variable bending 
rigidity at criticality lies in the Ising universality class (in the parameter regions II, 
III). From this we immediately infer that also the strong coupling lattice Schwinger 
model (rj = 1/2) belongs to this class. In || however a critical exponent v ~ 0.63 
has been reported for the strong coupling Schwinger model which is quite off the 
Ising result v = 1. The discrepancy very likely stems from finite size effects of the 
small lattices on non-quadratic domains investigated. These non-quadratic domain 
lattices however can be exploited in other ways as we will see in section 5. 



3 The independent loop approximation 

Inasmuch as exact expressions for the partition function (0) for general rj are not 
available analytical attempts to understand the phase structure of the self-avoiding 
loop model with variable bending rigidity have to rely on certain approximations. 
A method also applied in related situations in condensed matter physics and cosmic 
string physics is the so-called 'independent loop approximation' |l6l|-||2l|, p4| , |45 



The approximation is approached by writing the partition function Z A [z, rj\ as a sum 
over partition functions with a fixed number I of (polymer) loops. 



Z A [z,rj\ = z |A| £ Z A [l,z,rj\ (16) 

1=0 

The approximation now made is to express the Z-loop partition function Z A [l,z,rj\ 
exclusively by means of the single loop partition function Z A [1, z, rj] 

Za[1,z, V ] = i [Z A [l,z, V ] f , (17) 

leading to 



Z A [z, V } = z^ e^tl,*,*?] . (18) 

This approximation can be expected to give reasonable results for those parameter 
regions where the loop system is sufficiently dilute (filling in the average only a 
certain fraction of the lattice A). According to eq. ( |T8| ) the investigation now may 
concentrate on the single loop partition function Z A [1, z, r]]. One can easily convince 
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oneself that in the independent loop approximation the average number of loops in 
the system is given by the value of the single loop partition function |17]], eq. (56). 
The free energy density / reads in the independent loop approximation (f3 T = 1/T) 



(3 T f(z, V ) = - lim JL \nZ A [z, V ] = - In* - lim -(19) 

|A| — >oo | TV I |A| — >oo | TV | 

To proceed further, the single loop partition function can now be written as sum 
over the loop length 



Z A [l,z, V ] = £ z- 2k Z A [2k, V } (20) 

k=l 

where Z A [2k,r]] is the conformational partition function of a single loop of length 
2k (Here we already have taken into account that on a square lattice the length of 
a loop is always even.) pOfl . The conformational partition function is represented 
then as sum over all single loop configurations L with length 2k 



2k 

Z A [2k, V ] = £ rj c ^ = J2 N(2k,C) V c . (21) 

L, \L\=2k C=0 

Let us start with the consideration of the ordinary loop model (77 = 1). In this 
case Z A [2k, 1] denotes the total number of possibilities to place a self-avoiding loop 
with length 2k on the lattice A. It can be expressed by means of the number p2k 
of 2/c-step self-avoiding loops per lattice site which is a standard quantity that has 
been investigated in the literature. 



P2k 



lim 

|A| too 



Z A [2k,l] 
IAI 



(22) 



The n — > limit of the lattice 0(n) spin model provides us now just with the 
information necessary to study the critical behaviour |23[] (see also, e.g., 
2). For large k p 2 k reads |27| 



26 , sect. 



p 2k fe =°° B f k [2k]- 2 "- 1 + ... . (23) 

Here /1 denotes the connective constant (effective coordination number) for the self- 
avoiding walk problem on the given lattice A ^6|] and B is some lattice dependent 
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constant. The (universal) critical exponent v is believed to be given in two dimen- 
sions by v — | (obtained on a honeycomb lattice) f47j, j48|. Inserting ( p3|) into eq. 
one finds 



Z A [l,z,l] = \A\B J2 [2A:] 



-5/2 



fc=l 



(24) 



This is a justified approximation because we are mainly interested in the critical 
domain which is related to the k — > oo behaviour. From eq. (^4j) one easily recognizes 
that the critical point is given by z cr = fi. Most recent (precise) estimates for fi on 
the square lattice can be found in |]4~9|-|l51|. We keep here only a few digits and write 
z cr = A* = 2.638 (T cr = 1.031). Eq. (51) inserted into eq. (JT9[) gives immediately the 
free energy and one recognizes that the phase transition at z cr = \x = 2.638 found 
within the independent loop approximation is of second order. Using (fl52|, |53| ) 



F(x, k) 



OO n 

7jc 



n=l 



rr 



r(l - k) (-In a;) 



- 1 + E c(* 

n=0 



n) 



(In a;) 



(25) 



one reobtains for the critical exponent of the specific heat a the hyperscaling rela- 
tion a = 2 — 2u entailing in the independent loop approximation a = 1/2 which is 
to be confronted with the expected Ising result a 



CI. 



We are now prepared to study the general case with variable bending rigidity 
rj. First we have to find an appropriate expression for the number N(2k, C) of 
self-avoiding loops with length 2k and C corners. Let us count first the number 
of random non-backtracking walks of length 2k with C corners |2(J. It should be 
stressed that the following argument does not depend on the dimension of the lattice. 
There are 2k — 1 vertices available the C corners can be placed at, i.e., there are 
( 2fc (7 1 ) possibilities to do so. To each prospective corner exist h = (2d — 1) — 1 ways 
of bending where d is the dimension of a (hyper) cubic lattice (in our case of a square 
lattice d = 2). (2d — 1) is here the non-backtracking dimension of the lattice and it 
has to be diminished by 1 corresponding to the straight line choice. Consequently, 
we find 



N NB {2k,C) 



|A| 



'2k - r 

, c , 



(26) 



Using eq. ( ^1) the corresponding (non-backtracking) conformational partition func- 
tion reads then 
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Z A [2k, V ] NB = |A| [l + hr]} 2 "- 1 . (27) 

We now obtain an approximation to the 2/c-step self-avoiding walk (SAW) confor- 
mational partition function by simply replacing h = (2d — 1)[1 — l/(2d — 1)] by 
h = p[l — l/(2d — 1)] (this is based on the assumption that the self-avoidance con- 
straints effectively encoded in p are independent of whether propagation is straight 
or bent) where p is a certain effective coordination number to be determined in a 
moment. This yields in our case 



Z A [2k, V } 



SAW 



|A| U(2k) 



1 + 



2 m 



2fc-l 



(28) 



The additional factor U(2k) also to be determined below takes care of some ad- 
ditional length dependence which might show up in the transition from non-back- 
tracking to SAW. Specializing eq. ( p8|) to r\ = 1 we find the total number of SAW 
of length 2k on the square lattice which has to be confronted with the standard 
expectation pTfl for large k 



Z A [2k,l\sAW fe^oo a 2k roMT-l /nn\ 

c 2k = hm j-r-j = Aft 2/c 7 + ... (29) 

|A| — >oo | TV | 

where A is some lattice dependent constant and 7 = 43/32 ||7|], [fR|. From eq. ( |2"9"| ) 
we immediately find 



p = |(Ai-l) , (30) 

U(2k) = A p [2A;] 7 " 1 (31) 

Now, we need to know the conformational partition function for the self-avoiding 
loop problem. In order to be able to make further progress let us assume that 
N(2k, C) for self-avoiding loops is just a certain fraction M{2k) of N saw (2k, C) at 
least for large k independent of the number of corners C. According to eq. (|21|) then 
we can write 

Z A [2k, V } = M(2k) Z A [2k, V }sAw (32) 
which reads after having taken into account eqs. fl2"8]), fl3Ti|), ( |3"T| ) 
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Z A [2k,rj\ = |A| M(2k) A fi [ 1 + (/i - 1)77 



n 2fc-l 



[2*] 



7-1 



(33) 



M(2k) is the fraction to be determined. We here simply ignore the fact that for 
any loop the number of corners is even, necessarily. This is justified for the study 
of the k — > 00 behaviour we are primarily interested in. For rj = 1 we already 
have displayed an expression (eq. (p3"D) which now serves as reference expression to 
determine M(2k). We obtain 



M(2k) = j [2k}- 2u -~< 



leading to 



(34) 



Za[1,z,t]} 



|A| 



B /i 



[ 1 + (// - 1)77 j fc=1 
Consequently, the critical line is found to be 



E [2*] 



-5/2 



2 A- 



• (35) 



^7cr ( 2^cr ) 



(a*-i) 

This translates into the {T, s} coordinate system as 



(36) 



1/T a 



(37) 



It should be emphasized that the result of our approximate consideration (eq. (p7|)) 
entails s cr — > 1 for T cr — > 0. This is well in line with expectations spelled out in 
From the above equations we obtain for the strong coupling Schwinger model 
z cr = 1.819 (T cr = 0.774, s cr = 0.537). We find for the critical hopping parameter 



K cr (0) 



V / 2(^+l) 



0.371 



(38) 



The critical line (eqs. (p6|), (0)) obtained within the independent loop approxi- 
mation is plotted in figs. 4, 5. One recognizes that the critical line found analytically 
agrees qualitatively quite well with the result of the MC calculation of Miiser and 
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Rys pif . However, it is clear that the validity of the independent loop approxima- 
tion is confined to the low (polymer) loop density domain. The high density result 
of Miiser and Rys [^T|] displayed in region I of figs. 2, 3 cannot be obtained within 



the present scheme. It is also well known [16], [44] that while within the independent 



loop approximation the critical line can be determined in a qualitatively correct way 
results for critical exponents are correct to a lesser degree. This also applies to our 
case as we have seen above. While we would have expected, e.g., for the critical 
exponent a the Ising result (a = 0) within the independent loop approximation we 
see a = 1/2. Finally, it should also be mentioned that the relative simplicity of the 
independent loop approximation has its price because so far no way of its systematic 
improvement is known and one consequently has no quantitative control over the 
approximation made. Perhaps, this drawback is offset by the applicability of the 
approximation to systems in any dimension. 



4 The exact solution of the self-avoiding loop model 
with bending rigidity 77 = l/\/2 and its use 

While in general the self-avoiding loop model (SALM) with variable bending rigidity 
(which is equivalent to a 7-vertex model due to W2 = 0) cannot be studied exactly so 
far, there exists an exact solution to it for 77 = 1/Vz first investigated by Priezzhev 
[j37][ (see also f3~8fl ). This solution has later been rediscovered by Blum and Shapir 



who apparently were unaware of the earlier work of Priezzhev. The solution 



relies on the general study of the 8- vertex model by Fan and Wu p6[ . They 



found that the 8-vertex model is exactly solvable if the free fermion condition 



W1W2 + W3W4 



W 5 Wq + W 7 Wg 



(39) 



is fulfilled (cf. fig. 1 for the labelling of the vertices). Inserting eqs. ffl)-(0) into 
dUD (v taken arbitrary here) one immediately finds that the free fermion condition 
is fulfilled for 77 = 1/ y/2. The partition function for the SALM with bending rigidity 
r] = l/\/2 has been found in [[37]], by standard methods. The free energy density 
/ reads 



/W0M/V2 



I i-2-k (-2it 

' d6 dd> In 



8ir 2 Jo Jo 



2 + z 2 + 2,2 cos 6> 
+ 2z cos + 2 cos 9 cos <p } . (40) 
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A second order phase transition occurs for z cr = 2 (T cr = 2/(3 In 2) = 0.962, 
s CT = 1/3) which will be of main interest to us. There is of course also a criti- 
cal point at z = related to the exactly solvable 6- vertex model [52], ||. Because 



the system can be represented by means of free fermions [54] the SALM with bend- 



ing rigidity rj = 1/ v^2 lies in the Ising universality class . In accordance with 
this it has been shown (for w 2 chosen arbitrarily) that the partition function of the 



free fermion model can be expressed in terms of that of the regular Ising model |55 
It finally deserves mention that the result of the MC calculation of Miiser and Rys 
] j3lH is in complete agreement with the exact solution of the free fermion model (cf. 
figs. 2, 3). 

The above exact solution lying on the critical line of the SALM with variable 
bending rigidity is quite useful because this way one may take advantage of univer- 
sality arguments to draw conclusions about the model at criticality for fairly wide 
ranges of the bending rigidity rj. This will be discussed further in section 5. Here 
we will study the approximate calculation of the critical line in a neighbourhood of 
the model for rj = l/y/2. This discussion is in a certain sense a generalization of 



that given in [j38[ , f37| . Let us write the partition function (Q) in the following way. 



Z a [z,tj] = £ z^ 21 £ rfV (41) 

(=0 L, \L\=2l 



|A| 



1=0 L, \L\=2l k=0 



oo -y 



E^ |A| - 2 ' E £77[C(L)ln# (42) 



E E 

1=0 k=0 



E ^- 2l E W) a (43) 



Here, 



(C(L) k ) 2l = £ C(L) k . (44) 

L, \L\=2l 

One may now express {C(L))2i the following way. 

(C{L)) 2l = C 2l N(2l) (45) 
C 2l = 2ln c {2l) (46) 
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N(2l) = (1)21 denotes the number of (multi) loop configurations of total length 21 
on the lattice A and nc(2l), < nc{2l) < 1, stands for the average relative density 
of corners in the considered loop ensemble of total length 21. The following of course 
holds for the higher moments of C. 



< (C(L) k ) 2l < {2l) k N(2l) 



(47) 



One can now write 



Za[z,v] 



|A| 

E ^(20 
1=0 



"^n c (2Z)~ 


21 


Z 





({C{L)-C 2l f) 



21 



2 N(2l) 



(lanf 



+ 



(48) 



where ... stands for a series in higher order correlation functions of the corner 
number C and In 77. The critical behaviour of the system is related to large /. For 
I — > 00 nc tends to some value nc and consequently points { £1,771}, {22,^2} on the 
critical line not to far away from each other should obey to leading approximation 
the equation 



Tin 

Vi 



Z2 



(49) 



The contribution of correlation functions of C should be expected to be of minor 
importance in eq. (fi8"D, leading to corrections to the leading behaviour only. Inserting 
into eq. the exactly known critical point {z, 77} = {2, l/y/2} of the free fermion 
model leads in a neighbourhood of it to the equation for the critical line 



Vc 



-(n c +2)/2n c l/n c 



(50) 



The only unknown quantity in this expression is the average relative corner density 
nc- Its value nc(2l — > 00) is related to the high density polymer limit which is 

37[ 



reached for z — > 0. nc has been calculated in 



|38|j for the free fermion model 
and found to have the value 1/2. Values for the correlation functions of C for I — > 00 
can also be obtained along the same lines by tedious, but standard methods (fl5l 



|p7|1 ; the latter is the English original of ref. no. 16 in [37]). So, we end up with the 
following equation for the critical line of the SALM with variable bending rigidity 
in the neighbourhood of the free fermion point (cf. also figs. 4, 5). 
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Vcr {z cr ) = 2~ 5 / 2 z 2 cr . (51) 
This equation reads in {T, s} coordinates 

Consequently, we obtain for the strong coupling Schwinger model z cr = 2 3 / 4 « 1.682 
(T cr = 4/(7 In 2) « 0.824, s cr = 4/7 « 0.571). This yields for the critical hopping 
parameter 



/t cr (0) = 2~ 11/8 w 0.386 . (53) 

We see (cf. also figs. 4, 5) that the approximation based on the exactly solvable free 
fermion model yields a numerical value of the critical hopping parameter fairly close 
to the results of previous computer studies ]3"I] , ||, |J. As mentioned above, sys- 
tematic improvements can be obtained by taking into account correlation functions 
of C. This apparently is necessary as one learns from figs. 4, 5 if one wants to find 
the critical line beyond the region defined by the critical points of the ordinary loop 
model and the strong coupling Schwinger model respectively. 



5 Discussion and conclusions 

Let us first have a look at the larger picture emerging for the critical behaviour of the 
self-avoiding loop model with variable bending rigidity. There is one point on the 
critical line known exactly from the solution of the free fermion model (rj = l/v2, 



z cr = 2) ||37|-|]39l . For this model it is established that the phase transition is Ising- 
like, i.e., the model experiences a second order phase transition with exactly the 
same critical exponents as the regular Ising model. By the argument of universal- 
ity we may conclude that neighbouring models which lie on the same critical line 
exhibit the same behaviour. This in particular concerns the ordinary loop model 
(rj = 1) and the strong coupling Schwinger model (rj = 1/2). For the ordinary loop 
model this has been confirmed by MC investigations in the past f43| . For the strong 



coupling Schwinger model this consideration specifies the so far unknown character 
of the phase transition and confirms the recent suggestion that the transition might 
be a continuous one 101. 
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In order to extend the understanding of the self-avoiding loop model with variable 
bending rigidity at criticality let us consider the central charge c of the correspond- 
ing conformal field theory (CFT). Helpful information can be obtained most easily 
for the free fermion model considered in sect. 4. First, it seems worthwhile men- 
tioning that the regular Ising model can be understood as a special free fermion 



model [54]. A preliminary investigation along the lines given for the Ising model 
in [^] indicates that the self-avoiding loop model with bending rigidity 77 = 1/ v^2 
can be represented at the critical point z cr = 2 by just one massless (continuum) 
Majorana fermion (As in the special case of the Ising model just one half of the 
fermionic modes needed to express the partition function becomes massless at the 
critical point.). Consequently, this suggests that the critical self-avoiding loop model 
at the free fermion point is equivalent to a c = 1/2 CFT. The central charge cannot 
change continuously on the critical line in the neighbourhood of the free fermion 
model, therefore CFT's corresponding to the self-avoiding loop model with variable 
bending rigidity should all be expected to exhibit c = 1/2. This of course entails 
that the strong coupling Schwinger model at criticality should be equivalent to a 



c = 1/2 CFT. Consequently, in accordance with Zamolodchikov's c-theorem |59 



[ pCfl , z 7^ is related to a flow from the 6- vertex model (z = 0) having c = 1 Iplfl , 
towards a model with c = 1/2 (as discussed in general terms by Salmhofer M). 

We are going to test the above insight now by calculating the central charge for 
the strong coupling Schwinger model. This can be done most easily by considering 
the model on a strip of width a and length b — > 00 [p3H , [p4| . The central charge 
is related to the partition function (on a torus) by the formula (/ is the (bulk) free 
energy density on the infinite plane) 



= af( ^ + *1 (54) 

b^oo b 6 a 

We however will approach the study of the central charge of the strong coupling 
Schwinger model by means of the exact partition functions calculated on finite lat- 
tices in jH, ||. For our purpose the exact partition functions available on a 8 x 8 
torus 0], and on 2 x 32, 3 x 16 tori |5| are suited. In fig. 6 we have plotted the 
function 



c(z) = -- {\nZ A(axb) [z,l/2\ - abf(z, 1/2)} (55) 

for the 2 x 32, 3 x 16 tori in dependence on k (z = {k)~ 2 has been inserted). f(z, 1/2) 
has been calculated by means of the 8x8 partition function. For sufficiently large 
b the function c(z) should be expected to approach the value of the central charge 
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at the critical point. However, one has to be aware of the fact that on the very nar- 
row (with respect to a) tori considered massless and massive fields can contribute 
comparable amounts to the Casimir energy. Inasmuch as the central charge is cal- 
culated by means of eqs. (|54j), ([55]) from the Casimir energy results obtained from 
very narrow tori may turn out misleading. In part, this is what we observe from 
fig. 6. The result for the a = 2 torus (solid line) rather suggests c = 1 (or some 
value close to it), however the torus is so narrow that massless and massive fields 
contribute comparably to the Casimir energy. Consequently, in agreement with our 
expectation c = 1/2 for the wider a = 3 torus (dashed line) we already observe a 
much smaller value of c(z) at the critical point. However, it turns out that the sizes 
of the tori for which the exact partition functions have been calculated so far in the 
literature are too small to allow any final conclusions for the central charge of the 
strong coupling Schwinger model at criticality. Only partition functions calculated 
on considerably larger lattices will allow to numerically test the prediction c = 1/2 
in a reliable way. 

To conclude, the study of the self-avoiding loop model with variable bending 
rigidity presented in this article enhances the understanding of the critical behaviour 
of the strong coupling Schwinger model with Wilson fermions. We find that a second 
order phase transition which lies in the Ising model universality class takes place 
at some finite value of the hopping parameter K cr (0). Using certain approximate 
analytic methods the value of the critical hopping parameter is confirmed to lie at 
K cr ~ 0.38 — 0.39 in accordance with earlier numerical investigations [|J, ||. Certain 
arguments considered suggest that the strong coupling Schwinger model at critical- 
ity is equivalent to a c = 1/2 CFT. 

From a technical point of view, the present paper studies the application of the 
independent loop approximation to the qualitative and in part quantitative explo- 
ration of the phase structure of the self-avoiding loop model with variable bending 
rigidity in 2D. Comparison with known numerical results [|3T| shows that this method 
delivers a fairly correct picture for sufficiently low (polymer) loop densities. This is 
encouraging because the method is equally applicable to higher dimensions, while 
the analytic approach based on the exactly solvable free fermion model presented in 
section 4 is at least in part specific to 2D. This suggests that the independent loop 
approximation might successfully be applied also to analogous systems in higher di- 
mensions (e.g., to strong coupling QCD in 4D where the critical hopping parameter 
has recently been studied by other methods f6"5fl). 
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Figure captions 



Fig. 1 

Vertices of the 8-vertex model and their weights (cf. eqs. 
Fig. 2 

Critical line of the self-avoiding loop model with variable bending rigidity as found 
by MC calculations on a 64 x 64 lattice by Miiser and Rys |3l| . This figure has been 
copied and adapted from fig. 2 of |3"I| . The dashed line indicates some interpolated 
curve to the point s = 1, T = (specially emphasized by a black half disk). The do- 
mains I — VI are mapped to the correspondingly labelled domains in the {z, r/j-plane 
(see fig. 3, boundary lines are plotted in the same style in both figures), f f denotes 
the exactly known critical point {T cr = 2/(3 In 2) m 0.962, s cr = 1/3} (z cr = 2) of 
the free fermion model (77 = 1/V2) [p^]-|p9||. lg stands for the ordinary self-avoiding 
loop model (77 = 1, s = 0) with the critical point T cr = 1.157 [z CT = 2.373) j42jl-f43" 



sm denotes the critical point K cr (0) = 0.38 (T cr . = 0.80, s cr = 0.56, z cr = 1.73) of the 
strong coupling Schwinger model (77 = 1/2) as found in [|J, ||. 

Fig. 3 

This is the equivalent of fig. 2 shown here for the {z, 77} coordinate system. For 
further explanations refer to fig. 2. 

Fig. 4 

The critical line according to the results of the independent loop approximation 
( P7[) (dotted line) and of the free fermion model related approach (^) (dashed line) 
in comparison with the MC result (solid line) of Miiser and Rys |]3T|. For further 
explanations refer to fig. 2. 

Fig. 5 

This is the equivalent of fig. 4 shown here for the {z, 77} coordinate system and re- 
lating to eqs. ( ]3"6j ) (dotted line) and ( |5"T| ) (dashed line). 

Fig. 6 

The function c (see eq. (p5|)) in dependence on k(0). The solid line is the result for 
the 2 x 32 lattice while the dashed line stands for the result on the 3 x 16 lattice. 
The value of c at the critical point /t C r(0) — 0.38 — 0.39 has to be compared with 
the expectation for the central charge (For a discussion see the main text.) 
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